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Inspired by spin-electric couplings in molecular magnets, we introduce in the Kitaev honeycomb 
model a linear modification of the Ising interactions due to the presence of quantized cavity fields. 
This allows to control the properties of the low-energy toric code Hamiltonian, which can serve 
as a quantum memory, by tuning the physical parameters of the cavity modes, like frequencies, 
photon occupations, and coupling strengths. We study the properties of the model perturbatively 
by making use of the Schrieffer- Wolff transformation and show that, depending on the specific setup, 
the cavity modes can be useful in several ways. They allow to detect the presence of anyons through 
frequency shifts and to prolong the lifetime of the memory by enhancing the anyon excitation energy 
or mediating long-range anyon-anyon interactions with tunable sign. We consider both resonant and 
largely detuned cavity modes. 
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I. INTRODUCTION 

Realizing qubits with long coherence time is a ba- 
sic requirement for quantum information processing.^ A 
promising strategy to achieve this goal is to make use 
of quantum error correcting codes, where a single logical 
qubit is represented by suitable entangled states of a col- 
lection of many spins. This allows to detect and correct a 
finite number of errors acting on the individual physical 
spins of the codeP Furthermore, if the quantum code is 
the ground space of a gapped Hamiltonian, errors in the 
memory lead to excitations in the spin system and they 
can be suppressed at low temperature. In this respect, 
stabilizer Hamiltonians (given by a sum of mutual ly c om- 
muting Pauli operators) are of special importance^ and 
a particularly interesting class are stabilizer Hamiltoni- 
ans with topological order, where the distance of the code 
grows as a power of the system size N (see Refs. 0] and 
[5] for a rigorous definition of topological quantum order 
and its consequences on the stability of the ground sub- 
space under local Hamiltonian perturbations). Models 
with topological order wer e already known in the context 
of lattice gauge theory^] an d a pioneering example ap- 
plied to quantum information is the toric code proposed 
by KitaevP 

However, several limitations of this type of model 

appear when put in contact with an external 

environment PESl Most importantly, it was shown that 
at any finite temperature the lifetime o f the toric code 
does not grow with system sizePSE^EH This is a generic 
feature of stabilizer Hamiltonians with short-range inter- 
actions in one and two dimensions since, while the dis- 
tance of the code increases with N, the energy b arrier to 
perform a logical error is bounded by a constant! 15 * 16 ! In 
the case of the toric code, the excitations of the system 
can be represented as pairs of classical diffusing anyons 
and the aforementioned energy barrier is simply the en- 
ergy cost to create a single anyon pair. Therefore, it is 
of crucial importance to devise new architectures where 
such gap is enhanced and the anyon population is expo- 
nentially suppressed. 



It was recently shown that repulsive long-range inter- 
actions among the anyons lead to a large self-consistent 
mean-field gap which grows with system size and the 
resulting prolongation of the lifetime was studied in 
detailP^The toric code can be realized as a low-ener gy ef- 
fective Hamiltonian of the Kitaev honeycomb modelpES 
and long-range interactions between anyons appear in 
the presence of a non-local coupling with cavity modes 
extending over the whole memory!^ In analogy to the 
spin-electric coupling in molecular magnets ESHSS we con- 
sider a modification of the Ising couplings of the type 
Jx.y — > Jx,y + Sx,y(a> + a>). In this paper we study in 
detail the effect of such a coupling in the honeycomb 
model by making use of perturbation theory and show 
that cavity modes can be useful in several ways, not only 
to mediate long-range interactions between anyons. 

For example, they allow to realize a basic operation 
such as the read-out of the error syndrome or, in other 
words, to detect the presence of anyons across the mem- 
ory. To achieve this, it is sufficient to couple locally a 
single cavity mode to the desired read-out site and to 
detect its frequency shift, similar to a scheme demon- 
strated for superconducting qubitsP^Z] As we will see in 
the next sections, a single cavity mode can also be used to 
resonantly enhance the anyon gap while, to generate an 
anyon-anyon interaction of the desired sign, two cavity 
modes are required. Such a variety of situations corre- 
sponds to different specific configurations for the interac- 
tion between the cavity modes and the spin model, which 
we analyze by making use of perturbation theory. Our 
work is based on the Schrieffer- Wolff transformationpSHU 
which can be easily applied to the Kitaev honeycomb 
model, even in the presence of cavity modes, to obtain 
explicit expressions for the parameters of the quantum 
memory. 

The detailed outline of the paper is as follows. We 
start with Sec. |H]by introducing the model Hamiltonian 
and discuss the physical motivation of the coupling. We 
then describe in Sec. |III| t he perturbative framework in 
general terms. In Sec. |IV| we consider different read-out 
schemes. In Sec. [V] we study the enhancement of the 
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FIG. 1. Illustration of the Kitaev honeycomb model with 
cavity modes, a) Schematic view of the honeycomb model 
at the center of an optical or microwave cavity. The two 
shaded areas represent mirrors and the curves in the middle 
two quantized cavity modes with frequencies Wi,2 to which 
the spin model couples according to Eq. |2j. b) Kitaev hon- 
eycomb lattice. The spins are represented by white or black 
dots respectively belonging to the A or B sublattice. Ising 
interactions of type a x a x along x links (red), a y a y along y 
links (green), or a z a z along z links (blue) are represented by 
colored links with different directions. Strongly interacting 
z dimers are depicted by thick blue links and form a lattice 
with basis vectors ei,2 while t is a vector which connects A- 
and B-sublattice sites. Cavity modes (here, two of them, with 
frequencies wi.,2) modify the Ising interactions of certain links 
(as indicated by dashed lines), according to Eq. pj). 



gap produced by a single resonant cavity mode. In Sec. 
|VI| we obtain long-range anyon-anyon interactions from a 
specific coupling scheme with two resonant cavity modes. 
We conclude with Sec. VII| by studying the effect of off- 
resonant cavity modes and present in Appendices [A|P| 
an extended discussion of several technical aspects and 
derivations. 



II. MODEL 

The Kitaev honeycomb model 21 ^ is a two-dimensional 
version of the quant um compass model defined on a hon- 
eycomb latticePSED It consists of 27V spin-1/2 particles 



located at the sites of a honeycomb lattice with nearest- 
neighbor Ising interactions and is illustrated in Fig. [IJd) . 
The interaction between two neighboring spins is either 
a x a x , a v <J v , or a z a z , depending on the direction of the 
link connecting the two spins. In this work we are in- 
terested in the limit when the a z a z couplings are much 
stronger than the other ones. Therefore, the model can 
be seen as consisting of N weakly coupled dimers, form- 
ing a triangular lattice with unit vectors ei 2 , as shown 
in Fig.[T],;.. 

In this work, we consider the following generalization 
of the Kitaev honeycomb Hamiltonian 



H = 



E 



K^K<+ ei +t - ^K^ y alal +e2+t , (1) 



where the three last sums run over all vectors a = 
n\e\ + (n\,2 G Z) pointing to a ^4-sublattice site 

(white dots), t is a vector which connects A- (white dots) 
and B- (black dots) sublattice sites and a Xl y ' z are the 
usual spin Pauli operators. In Eq. ([I]), we have assumed 
that the Ising couplings K a j, (with k — x,y, z) are site 
dependent. Furthermore, we introduced a collection of 
cavity modes labeled by a, with frequency u) a , and anni- 
hilation operators a a . We propose that they are coupled 
to the spin model as follows: 



(2) 



where J a ^ are the unperturbed Ising couplings, and 
S a ,k,a describe the linear change due to the cavity modes 
(see Fig. [I). The original Kitaev model is recovered for 
<5a,fc,a = an d J a j. = Jfc, independent of the lattice site. 
For simplicity, we also assume in the rest of this work 
K z = J z , but we keep the more general form of coupling 
[Eq. ([1] for K a , x and K a , y . 

To physically justify Eq. ([2|, we note that Ising cou- 
plings can be engineered in a variety of systems, for 
example between quantum dotsp^l Josephson Junction 
qubitsj22H2U ; atoms in optical latticesJ^S^ZI anc l doped 
coupled cavities. 38 Anisotropic spin interactions are also 
present in molecular magnets^HlSl and between pseu- 
dospin orbital states.^ Recently, theoretical proposals 
to engineer the honeycomb model in Mott insulators 
with strong spin-orbit coupling were discussed in Refs. 
I5T1 and Metal-oxide compounds (layered iridates of 
type J42RO3, with A = Li,Na) were considered there 
as promising materials. In these systems, the strength 
of the magnetic interaction can be modified by exter- 
nal per turbatio ns, in particular electric and magnetic 
fieldsPt25132l34] g uch perturbations can be generated by 
charged nanomechanical systems (e.g., cantilevers) or 
quantized electromagnetic fields (e.g., in cavities and 
transmission line resonator d 26 * 27 * 40 *) , thus realizing a cou- 
pling of the type of Eq. As a more specific example, 
we would like to mention molecular magnets. In this 
case, electric fields can be used to modify the overlap of 



FIG. 2. Pictorial representation of superexchange interac- 
tions. The bridge sites (dots in the middle of each i-link) 
mediate the magnetic a x a x interactions between the spins 
(dots at each site of the honeycomb lattice). Here, the x links 
have a nonvanishing dipole moment and can thus be linearly 
coupled to external electric fields. By controlling the equilib- 
rium position of the bridge sites, the coupling to the cavity 
modes can in principle be tuned locally. 



the orbital wave functions and thus the resulting mag- 
netic interactions, which are determined by exchange or 
superexchange mechanisms. Strongly anisotropic (Ising- 
like) magnetic interactions can arise here in the presence 
of spin-orbit coupling pHHU 

A linear effect in the electric field, as in Eq. j2j, can 
only exist if the electric-dipole of the underlying bond is 
not zero or, in other terms, if the inversion symmetry of 
that bond is broken. This is illustrated in Fig. [2] for the 
x links, assuming the presence of superexchange interac- 
tions mediated via an auxiliary site (the so-called bridge 
site) . Figure [2] also illustrates that only selected links 
can be coupled to the cavity modes since y- and z-bonds 
have no dipole moment and thus <5 a ,fc,a = 0. Finally, 
Fig. [2] suggests that 5 & ^,a can be controlled locally. By 
realizing Fig. [2] in a quantum dot setup, 5 a ,fe,a could be 
made zero or not by changing the equilibrium position of 
the bridge sites with static electric gates. 



III. PERTURBATIVE APPROACH 

Exact results for the Kitaev honeycomb model exist for 
arbitrary values of J x , y ,z 21 but they do not include cou- 
plings to the cavity modes. Therefore, we study Eq. ([I]) 
with perturbation theory. This is indeed appropriate if 
the model is intended to be a quantum memory, since 
for J x y <C J z the h oneycomb model can be mapped to 
the toric codel^HSU Without cavity modes, perturbative 
treatments of the honeycomb model can be found in Ref. 
[21] where the author makes use of Green's functions, and 
more recently in Ref. [22] with the p erturb ative contin- 
uous unitary transformation method! 41 ! 42 ! This work is 
based on a different approach, namely on the Schrieffer- 



FIG. 3. Mapping of the Kitaev honeycomb model to the toric 
code. Left-hand side: each z dimer, formed by a black and 
a white spin [see also Fig. [lj>)] , corresponds to a new lattice 
site (large blue dots). To these new sites, hard-core bosons 
and effective spins are associated. A low-energy (zero hard- 
core bosons) toric code Hamiltonian for the effective spins 
is obtained by perturbation theory. The resulting toric code 
model is schematically illustrated on the right-hand side. The 
dark and white unit cells on the left-hand-side correspond now 
to plaquette (p) and star (s) operators, respectively. 

Wolff (SW) transformationplMI which can be applied in 
a straightforward way to Eq. ([I]), as we will see in the 
next sections. 



A. Hard-core boson transformation 

The first step of our analysis consists in the hard-core 
boson transformation presented in Ref. 1221 which maps 
each z dimer to a particle and an effective spin, see Fig. [3] 
The two spins of each z-dimer can be either parallel or 
antiparallel and, if we assume that J z > 0, the parallel 
configuration has lower energy than the antiparallel one. 
The excited dimers can thus be in terpreted as particles 
(hard-core bosons) with energy 2J Z , and the remaining 
degree of freedom as an effective spin- 1/2. More explic- 
itly, 

I t)a ® | t)a+t, = I it 0)a, U)a ® | l)a+t = | 4 0) aj 

I Da® | t)a+t = I t l)a, | t)a ® | l)a+t, = | ^ l)a, 

where the four possible spin configurations of a given 
dimer (left-hand sides) are mapped to states with 0, 1 
hard-core boson and 1|, fr effective spin. If we rewrite the 
above transformation in operator language we obtain 

< + t=T*(bl + b a ), < = & a + 6a) 

°l+t = *M + b a ), at = iT*M - 6 a ), (3) 

<+t=r*, *S = t£(1 - 26£& a ), 

where 6 a (& a ) creates (destroys) a hard-core boson and 
T x,y,z are p au ji operators of the effective spin-1/2 at 
site a. At distinct sites a and a', bosonic commutation 
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relations are satisfied: 



(4) 



Furthermore, {& a -fr a } = 1. 

By making use of these bosonic operators, we can 
rewrite the first line of Eq. ([!]) simply as 



(5) 



notation introduced in Eqs. (7|9!, which define the f n 
functions, we can write 



— /n(Ja,fc)) 
^ Q — fn{5 a .k 



fn (^a,fc,a)&a 



(14) 
(15) 
(16) 



Clearly, the T n is also given by a sum of the T™ a : 



which represents the unperturbed Hamiltonian. Here and 
in the rest of the paper we drop an inessential constant 



B. Perturbation operators 

We now consider the second line of Eq. ( lj , which con- 
stitutes our perturbation. By rewriting the spin opera- 
tors by making use of Eq. (J3|, we obtain that Eq. ([I]) 
takes the following form 



H = Ho + T Q + T- 2 + T. 



+2, 



where 



T = - J2( K *M*,k + H.c.) := f (K, 

EL.k 

T-2 = {T +2 ) ] := f-2(K a ,k), 



(6) 

(7) 
(8) 
(9) 



with k = x,y. As in Ref. we defined the hopping 
operators 



-ibl 



h t x 
-ei u a'a+ei ' 

-e 2 "a ' a +e 2 'a > 



and the hard-core boson creation operators 

v — frl T x 

u a.,x — u a+ei"a'a+eiJ 

Va, y = ib\^„birl_ 



a+e 2 a a+e 2 a ■ 



(10) 

(11) 

(12) 
(13) 



In other words, T +2 (-2) creates (destroys) two nearest- 
neighbor hard-core bosons and To makes a hard-core bo- 
son hop to a nearest-neighbor unoccupied dimer. 

It is also useful to keep track of the change in the num- 
ber of photons: whenever a hard-core boson is created, 
destroyed, or hops from one site to another, a photon 
might be simultaneously created or destroyed. There- 
fore, we introduce T™ a operators, where the lower label 
(n = 0, ±2, as before) refers to the change in the number 
of hard-core bosons, while the upper label (m = 0,±1) 
indicates that the number of photons in mode a changes 
by m. Such operators are simply defined by substituting 
in Eqs. (7|9) the full couplings K a k by the appropri- 
ate quantity, J a<k , 5 & ^ a a a , or (5 a ,fc, Q aJ,. By using the 



T n = T° + Y,Ti 



(17) 



a,± 



Furthermore, the energy change corresponding to T™ a 
is immediately obtained from the energies 2 J z and uj a of 
the hard-core bosons and the photons respectively, and 
can be expressed through the following commutation re- 
lations 



[H^T™ a ] = {2nJ z +mw a )TZ 



(18) 



Note that we use here the convention 0a = such that 
Tg and 7± 2 , which do not change the state of the cavity, 
are defined independently of the value of a. 



SW transformation and correspondence to the 
toric code 



In the perturbative limit we are interested in, the SW 
transformation allows to obtain an effective Hamiltonian 
in a desired subspace up to an arbitrary order in per- 
turbation theory. For the convenience of the reader, we 
summarize the general procedure in Appendix |Aj and 
provide there the general formulas appropriate for the 
type of Hamiltonian of interest in this work, up to the 
fourth perturbative order. Since the quantum informa- 
tion is encoded in the spin degrees of freedom, we always 
consider the low-energy subspace where no hard-core bo- 
son is present. On the other hand, we will generally allow 
the modes of the cavity to be excited, to study how the 
presence of a finite photon population affects the proper- 
ties of the memory. While the treatment of Appendix [S] 
is completely general, the resulting fourth-order expres- 
sions involve too many terms to be presented here. It is 
instead interesting to consider specific coupling schemes 
and physically relevant regimes, in which only a few dom- 
inant contributions are important. Several examples will 
be examined in detail in the following sections. 

In the remainder of this section we restrict ourselves to 
the case without cavity modes, and derive the toric code 
from the formulas of Appendix|Aj By setting all the <J a ,fc,c« 
to zero, the T n operators coincide with T° and all the 



summations on the photon indexes i, j, k, r in Eqs. ( A24 
A26 1 can be dropped. By applying such formulas in the 
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FIG 

K, 



4. Honeycomb lattice with two distinct ^-couplings 
X1 2 and two distinct j/-couplings K yi 2 . 



and can thus be dropped for the moment. We have kept 
them here since they will become important in Sec. VII 
when the J X12 , J Vl 2 will be generalized to full couplings 
K Xl 2 , K yi 2 , including the cavity modes. 

The last line of Eq. p0| ) is expressed in terms of the 
spin operators W a , defined as 



W a 



T z y y z 

a a+e2 a+ei a 



(21) 



Note that the W a operator involves four dimers, con- 
nected by two x links and two y links. We will use 
sometimes the expression that such links "belong to" the 
corresponding operator W a . The W a operators are mu- 
tually commuting and can be rewritten in terms of stars 
and plaquettes of the toric code, as introduced in Ref. |H1 
In order to see this correspondence, we perform a spin 
rotation.^ If the spin lies at a site which is the bottom 
corner of a dark unit cell (see the left-hand side of Fig. [3]) 

t x = -s y , r y = s x , t z = s z , (22) 



subspace with zero hard-core bosons we obtain 



ir z T - 2T+2 + vhi T ~ 2ToT+2 

rri rri m rri rri rrt rri rri 

J--2-L-2-L+2-L+2 ~ n , — 2 ~\~2 



128JJ 
1 



64 J? 



64JJ 



T-2T+2T-2T+2, 



(19) 



where we wrote explicitly the third-order term in the 
second line, even if it gives no contribution: a hard-core- 
boson pair created by T +2 from the vacuum, followed 
by a To hopping process, cannot be annihilated by the 
T_2 operator. In general, only even orders appear in the 
perturbation theory and Eq. (A25l always evaluates to 
zero. It is also worth pointing out that Eq. (19) is the 
same as the one derived in Ref. [55] with the perturbative 
continuous unitary transformation approach. 

If we now explicitly write Eq. ( 19 1 in terms of spin 
operators, by using Eqs. ([7][9| with K ak = J k , the toric 
code Hamiltonain is obtaineaP As a slight generalization, 



we consider here the case when four distinct couplings 
enter, namely J Xl 2 and J yi 2 . This specific scheme is 
illustrated in Fig. [4] and leads to 



H r fr — 



I 2 



J x 1 J y 1 



4J Z 
j2 j2 



~ 2 
N 



J yi 



Jy,N 



AJ Z 

j2 j2 



r r 

J x 2 J yi 



N 



16J 3 



p p 

U X! °X2 



j2 p 

J yi J y2 



2 

N 



J X! 



16J| 



16J3 



64J 3 



+ 



64J3 



16J 3 



(20) 



This result will be useful later, when in Sec. |VII| we will 
discuss the non-resonant coupling of the cavity modes to 
the spin Hamiltonian. The first three lines are a constant 



otherwise 



(23) 



We have now that half of the unit cells of the lattice 
of dimers (the dark unit cells in the left-hand side of 
Fig. [3]) correspond to plaquette operators B p of the form 
s z s z s z s z . The other half of the cells (the empty ones in 
the left-hand side of Fig. [3| correspond to star operators 
A s of the form s x s x s x s x . By setting J Xl = J X2 = J x and 

and i2J 



J y , Eq. (20) gives back the result of Refs. [5] 



H. 



off 




(24) 



In the following we will use the simpler notation in terms 
of the W a operators, by keeping in mind that 

Y,w a = Y,A s +Y. R p- ( 25 ) 



IV. READ-OUT SCHEMES 

As a first application of the coupling of the honeycomb 
spin model to cavity modes, we show that it allows the 
read-out of the eigenvalues of the star and plaquette oper- 
ators. This is a basic requirement if the toric code model 
is aimed to be a quantum memory, since the knowledge of 
such eigenvalues (the error syndrome) is needed to per- 
form the error-correction algorithm and to retrieve the 
encoded information. 

We start by discussing an example when a single cav- 
ity mode, with unperturbed frequency u) and annihilation 
operator a, allows to read the eigenvalue of a single op- 
erator Wa r , for a given site a r . The specific coupling 



G 



scheme we consider for Eq. ^ is as follows, for the 
links: 



a. r +e 2 ,x 



= Jx +S x (a + a f ), 



(26) 



and K &tX — J x otherwise. Furthermore, K a y = J y for 
all the y links. In other terms, only the two x links 
which belong to the W &r operator are coupled to the 
cavity mode. This allows to read out the eigenvalue of 
W &r from an indirect measurement since we obtain the 
following effective Hamiltonian: 



H, 



251 



46U? 



4J Z - UJ (4J Z - uj} 



era, 



(27) 



which is derived from Eq. (JlJ by applying the general 
formalism of Appendix [A} In particular, as described in 
more detail in AppendixlB] we have evaluated Eqs. ( A24) 



and ( A26 ) for the specific coupling scheme under consid- 



eration, keeping only the terms diagonal in the number 
of photons and assuming the resonant condition 



(4J Z - w) < 4J z ,w. 



(28) 



This implies that, at each order, only resonant terms con- 
tribute to Eq. (27), i.e., the ones where only powers of 
(4 J z — uj) appear in the energy denominators. Here, and 
in the rest of this work, we also assume (4J Z — a;) > 0, 
such that the state with zero hard-core bosons is the un- 
perturbed ground state of the resonant subspace. 

Equation (27) shows that the frequency of the cavity 
is now 



251 



UJ 



4 J z - oj (4 J z - uj) 3 



(29) 



Therefore, it is sufficient to detect the change in the fre- 
quency of the cavity mode, which is conditioned by the 



eigenvalue of W &T {±1). Let us also recall from Sec. IIIC 
that the lattice site a r determines if the operator W &r , 
defined in Eq. plj ), corresponds to a star (empty unit 
cells in the left hand side of Fig. [3]) or a plaquette (dark 
unit cells in the left hand side of Fig. [3| of the toric code 
model. 

It is worth mentioning that a similar read-out scheme 
was demonstrated in superconducting circuits for a sin- 
gle qubit: in that case the role of W &r is played by cr z of 
the qubit being re ad, di spersively coupled to a transmis- 
sion line resonator! 26 * 27 ! The frequency shift of the cavity 
mode can be determined in transmission measurements. 
Indeed, proposals for the realization of the honeycomb 
spin Hamiltonian with superconducting circuits already 
exist pSl and it is conceivable that our readout scheme 
could also be realized in this type of systems. Note that 
the read-out must take place in a time interval which is 
much smaller than the typical time for a spin to flip. This 
condition can be satisfied for a single superconducting 
qubit, since the relaxation time can reach about 10 [is, 
while the response time of the measurement is a few hun- 
dred nanoseconds.^ However, the signal to noise ratio 



is presently too low to allow single-shot read-out. Be- 
sides superconducting systems, also spin qubits in semi- 
conductor nanostructures can be coupled to transmission 
line resonators via spin-orbit interactions. In particular, 
InAs nanowire quantum dots were discussed in Ref. 1431 

We now address the problem of measuring all the stars 
and plaquettes of the memory. As discussed already at 
the end of Sec. |Hj we can imagine that all the couplings 
5 a ,x,a can be switched on/off. For our read-out scheme, 
this implies that the site a r can be freely chosen and all 
plaquettes and stars can be successively measured, one 
after the other. This method is not efficient in terms of 
the total read-out time for the whole memory, but has the 
advantage of using a single cavity mode. Having at dis- 
posal many cavity modes (ideally, one for each site of the 
memory), much more efficient schemes can be realized. 

Consider in particular a set of sites £1 = {a Q }, each 
of them corresponding to an operator W aa which we 
would like to measure with a cavity mode a. For sim- 
plicity, we assume that all frequencies and couplings are 
equal, which can be realized with an array of identical 
resonators. Each mode a (in particular, the read-out of 
its frequency) can be addressed individually and is cou- 
pled to the two x links of the corresponding W &ol operator 
in the following way 

K &a , x = K aa+e2 , x = J x + 5 x (a a + a'a), (30) 

while K a x — J x and K a y = J y for all the remaining 
links. We further require that 



if a Q € fi, then a Q ± e2 ^ ft, 



(31) 



such that two operators W a , W &a2 never share a x-link 
and thus no virtual process couples the W &tx operators to 
more than one cavity mode in lowest order. Then, the 
read-out of a single site can be generalized to obtain the 
following effective Hamiltonian: 



251 



AS*J* 



4J Z - uj (4J Z - uj) : 



:W a 



(32) 

Therefore, we can read-out simultaneously the eigenval- 
ues of all W a „ through the shifts of the corresponding 
cavity mode frequencies. This procedure can be used to 
measure the state of the full memory in two steps only. 
This requires to have a cavity mode for each W a op- 
erator and being able to turn on/off all the spin-electric 
couplings. Then, we can apply the above procedure twice 
in a checkerboard configuration where, in the first step, 
we turn on the couplings of only the plaquette operators 
and, in the second step, the couplings of only the star 
operators. 



V. RESONANT ENHANCEMENT OF THE GAP 
FROM A SINGLE CAVITY MODE 

We discuss in this and the following sections how the 
cavity modes can be used to prolong the lifetime of the 
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encoded information. As shown in Sec. |III| the honey- 
comb model ([!]) in the absence of cavity modes and with 
uniform couplings (K ax = J x and K M = J y ) is equiv- 
alent to the toric code model of Eq. (24), up to fourth 



order in perturbation theory. This effective model can be 
used as a quantum memory encoding two qubits in the 
ground space, and the lifetime is essentially de term ined 
by the energy gapP As it can be seen in Eq. (24), the 
excitations correspond to sites where W a = —1. They 
can also be interpreted as anyons,^ with anyon numbers 
n„ defined as 



W a = 1 - 2n„ 



at each site a of the lattice. If W a = +1, then no anyon 
is present, i.e., n a = 0. Instead, if W a = — 1 one anyon 
is present, i.e. n a = +1. In Eq. (24), the anyons are 
nonintcracting particles with energy 



Mo = 



72 t2 

°x J y 

8J? ■ 



(34) 



In the present section, we consider how to obtain a non- 
interacting toric model in the presence of a single cavity 
mode, and how the gap of the model is affected by the 
cavity. We refer to the next section for a discussion of 
anyon interactions induced by two cavity modes, as pro- 
posed in Ref. [TBI 



As for the single-site read-out discussed in Sec. IV the 
cavity mode has unperturbed frequency w and annihila- 
tion operator a, but the coupling scheme is as follows, for 
all the x links: 



K a , x = J x +5 x (a + a f ), 



(35) 



while K„ 



J y for all the y links. Our analysis can 



also be simply extended to non-homogeneous couplings, 
when 5 X becomes site dependent. The read-out scheme 
of Eq. (26), where only two x links are affected by the 



cavity, provides an explicit example. 

To obtain an effective Hamiltonian, we assume again 
here the resonant condition of Eq. ( 28 1 , which allows us to 



keep only the leading term (with resonant energy denom- 
inators) in the perturbative contributions of Eqs. (A24) 
and (A26). The interested reader can find more details 



of the derivation in Appendix ml The final result reads 



51 N 



off 



4J Z 



-a) a- 



(36) 



where the modified gap of the model is immediately seen 
to be 



M = 



{AJ Z - u)f 



(a f o), 



(37) 



which is proportional to the number of photons populat- 
ing the cavity mode. In fact, comparing \i to the bare gap 
/j, of Eq. ( 34 ) , the quantity d x (a 1 * a) appears instead of J x . 
Furthermore, the denominator is in terms of (4J 2 — u>), 



which can be made smaller by simply changing the fre- 
quency of the cavity, in contrast to the bare value J z . 
Indeed, we assumed (4J Z — uj) <C J z , see Eq. (28). These 



facts make the effect of that cavity very interesting, be- 
cause the resonant gap fi can be made larger than hq. 
Here and in the rest of this work, we assume that the 
cavity is driven out of equilibrium. The average (a^a) 



appearing in Eq. ( 37 1 corresponds to the nonequilibrium 



photon population of the cavity and can be large. 

Note however that the gap cannot be made arbitrarily 
large. In particular, the perturbative expansion requires 
5 X ^/ (at a), J x y < (4J 2 — w), which implies /i < J y . Still, 
(33) for given values J x _ y <C J z , an enhancement factor of 



the bare gap /iq of order J z /J x J y can be achieved, by 
an appropriate design of the cavity and excitation of the 
resonant mode. Such increase in the gap has a dramatic 
effect on the lifetime of the quantum memory, if fi can 
become larger than the temperature T. In this case, the 
population of anyons is exponentially suppressed and the 
lifetime r increases accordingly. The following approxi- 



(38) 



mate formula was propose d 12 * 14 * 18 !: 



2/c 



e ji/k B T , 1 



D 



where D is the diffusion constant of the anyons and the 
prefactor f c can be interpreted as a critical fraction of 
errors.^ The value of D depends on the details of the 
thermal bath, especially whether it is Ohmic or super- 
Ohmic, and can also contain a dependence on fi (see Ref. 
[TH|for an extended discussion). 

Finally, we recall that the exact solution in the ab- 
sence of cavity modes^ also displays a gapless phase, 
away from the perturbative regime (i.e., when J x ,y,z have 
comparable values). The influence of the cavity mode on 
this gapless phase is an interesting question: Since in the 
resonant case the role of 4J Z is played by 4J Z — ui, it can 
be expected that not only the gap of the system can be 
modified by tuning u) around 4J Z , but also a transition 
from the gapped to the gapless phase could be induced by 
changing the parameters of the cavity. These nonpertur- 
bative aspects will be the subject of future investigations. 



VI. LONG-RANGE INTERACTIONS 

It was shown in Ref. [18] that repulsive long-range in- 
teractions between anyons have a beneficial effect on the 
memory. Since cavity modes spatially uniform across the 
memory realize interactions with constant strength, we 
consider the following effective model: 



H ci f = M X/' 



.4 



\ ^ 



2 



(39) 



where a, a' run over all the N sites of the lattice^ and 
the anyon numbers n a = 0, 1 are defined as in Eq. p3). 



Note that, differently from Ref. [JJ3 we do not distinguish 
here between star and plaquette anyons, but all of them 



interact among each other. The effect of the interactions 
on the memory lifetime can be understood in terms of a 
mean- field gap AW, which includes the repulsion energy 
(if A > 0) from the average anyon density n m f: 



AW = M + n ml AN, 



(40) 



where n m { is determined self-consistently as n m { — 
l/(e Atmf / fcBT + 1) assuming an equilibrium density of 
anyons at temperature T. 

It is the purpose of this section to study how the model 



( 39 1 is realized perturbatively and when the effect of the 
interactions becomes relevant. A noticeable feature of 



Eq. (40) is that, assuming constant fi, A, the mean-field 
gap increases with N. This is because the contribution 
from the interaction (second term) grows and becomes 
eventually the dominant term for sufficiently large N. 
The effect is to reduce the anyon density n m f and pro- 
long the memory lifetime.^ However, two main differ- 
ences appear in the explicit perturbative derivation: 

(i) The parameters \i and A acquire in general a non- 
trivial dependence on N. Therefore, it is possible that 
the anyon gap /i dominates the size dependence of the 
gap, instead of the interaction contribution. 

(ii) The requirement to strictly remain in the perturba- 
tive regime imposes restrictions to the size of the system, 
below which in many cases the interaction contribution 
is small with respect to /i. 

However, we also find specific coupling schemes and ap- 
propriate range of parameters for which the interactions 
become the dominant effect. 



A. Coupling scheme and resonant effective 
Hamiltonian 

As suggested in Ref. [TSl we consider two cavity modes 
with frequencies ui x ,y and annihilation operators a x _ y . To 
have two modes is useful because it allows to change the 
sign of the interaction, by choosing the frequencies and 
photon occupations of the two cavity modes. The first 
(second) mode couples to one-half of the x(y)-couplings. 
In particular, by referring to Fig. |1J we assume 

K^x = Jx+S x (a x +a x ), for the dark (red) x links, (41) 

while K ax = J x for the light (yellow) x links and 

K a y = Jy + 6 y (al + CLy) , for the dark (green) y links, 

(42) 

while K ay = J y for the light (white) y links. 

We obtain the effective Hamiltonian for the toric code 
model from a SW transformation by assuming, as in the 
previous sections, that the two modes are resonant: 



(4J Z 



(43) 



Although this considerably simplifies the treatment, still 
several contributions are present when evaluating the sec- 
ond and fourth order expressions, as discussed in more 



detail in Appendix [C| The following much simpler ex- 
pression is obtained: 



-ffeff = u' x ala x + J y a\a y 



J xJy &x ^y 
^y i^-^) 3 



E 



(4 a y + a x a y ) ^ W B , (44) 



by imposing the more restrictive condition 

\uj' x - u' \ < (4J Z - u) xy ) < J z ,u) x> . 



In Eqs. (44 1 and (45) we defined 



N 



2 4J Z 



(45) 



(46) 



where the second term is generally a small correction to 
the frequencies, due to perturbative restrictions to the 
size of the system [see Eq. ( B8 1] . A feature of the scheme 
considered here is that no term diagonal in the photon 
modes which couples to ^ a W a appears at fourth or- 
der. This is because each cavity mode only interacts 
with a single link of each W a operator, and thus at most 
one photon operator per mode can appear in combina- 
tion with ^ a W a . Therefore, all the couplings to W a in 
the effective fourth-order Hamiltonian are off-diagonal in 
the photon operators. As a first approximation, due to 
Eq. (45 ), we only kept one of these terms, namely the one 



appearing in the second line of Eq. ( 44 ) which is resonant 



in the difference of the two cavity mode frequencies. 



By performing a second SW transformation of Eq. ( 44 ) 



in the photon modes (in the same way described in Ref. 
I18p and keeping only terms involving spin operators, we 
obtain 



where A is is given by 



A 



A = 2 



E 



JxJy $x $y 



]_ct=x,y v z 



W Q ) E 



w a 



(4a x ) - {oydy) 



(47) 



(48) 



By rewriting Eq. ( 47 ) in terms of anyons we obtain 

AN ^ A 

H e $ = — n a + — n a n al . 



(49) 



As announced, the effective parameters /x and A have 
a dependence on TV which in the case of the interaction 
strength A can be quite weak [it a ppears through the 
denominator uj' x — uj' see Eqs. (48) and (46)]. Instead, 



the chemical potential is approximately linear in N 



A* = 



AN 



(50) 



which makes it the dominant effect. Therefore, it is re- 
quired for the stability of the n a = ground state that 
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A < 0, which leads to /i > and can be realized by an 
appropriate choice of the frequencies oj' a and photon oc- 
cupations (a^cia). The anyon interaction is in this case 
negative, which has the unfavorable effect of reducing 
the noninteracting gap, but it becomes quickly negligible 
with TV. Note that, for a gap which grows linearly with 
system size, the lifetime is prolonged exponentially with 
N at low temperature, as shown in Eq. (38). 



C. Off-resonant contributions to the anyon gap 



We have assumed so far the resonant condition ( 45 1 , 



such that we could keep only a single term which cou- 
ples to the star and plaquette operators [the second line 
of Eq. ( 44 )] . However, many contributions to the gap 
and interactions exist in general and we show here that 
for some regime of parameters it is possible that (i in 
Eq. (39 1 becomes zero. In this case, the interactions be- 



B. Interpretation in terms of anyon holes 



tween anyons have a dominant effect. 

We consider again the effective photon Hamiltonian up 
to fourth order, obtained from Eqs. ( A23) and (A26), but 



An interesting aspect of Eq. (|47) is that it is symmetric keep two additional terms with respect to Eq. 1441 



with respect to the change of sign of all the W a . There- 
fore, it is useful to define the anyon-hole numbers 



1 



(51) 



which describe excitations of the ground state with W a = 
— 1 for every a. Clearly, transforming Eq. (47) in terms 
of the n a leads to an Hamiltonian with the same form of 
Eq. ( 49 1 . Another interesting way to rewrite Eq. ( 47 1 in 



a symmetric way is by considering both types of particles 
present in the memory (i.e., at each site either an anyon 
or an anyon hole is present). This gives 



(52) 



describing a long-range interaction between anyons and 
anyon holes. For A < such interaction is repulsive and 
the ground state is completely occupied with one type 
of particles, say anyon holes. The appearance of anyons 
in the memory then results in a repulsive contribution 
from the ~ N anyon holes already present in the ground 
state. This picture provides a natural interpretation of 
the system size dependence of the gap, in terms of long- 
range interactions. 

As a side remark, we consider again the model Hamil- 
tonian Eq. (39), but with with constant coefficients \i,A. 



By rewriting Eq. ( 39 ) in terms of anyon holes we obtain 



H eS = -(n + AN)'£ i >, 



A^_ _ 

a. a' 



(53) 



Differently from Ref. Q21 we now assume A < 0. Then, 
the first term of Hamiltonian ( 53 ) is a size-dependent 



chemical potential which for large system size becomes 
positive, while the second term describes a long-range 
attractive interaction. Therefore, if we want to use this 
system as a quantum memory, we can encode a state in 
the memory full of anyons (W a = — 1 for every a, instead 
of +1). In this situation, the role of the anyons and the 
anyon holes is interchanged. For example, a spin-flip will 
now produce a pair of anyon holes which diffuse in the 
memory and destroy the information. However, since the 
chemical potential of the anyon holes increases linearly 
with the system size, one can prolong the lifetime of the 
memory by making the system size bigger. 



H eff = u; x a x a x + u y a^a y - - ^ AT _,, ( a « + a ") 

a—x,y 



J 2 T A 

E' J c\ ' J a u a 
UJ Z -w Q ) 3 



(a a + a*) 



E^ 



E 



2{U z ~u: a f 



(ala y + a x al)J2w a , (54) 



where in the second line a = y if a = x and vice versa. 
For some details on the derivation of this expression, we 



refer to Appendix \Cl In Eq. ( 54 ) , the second line gener- 



ates now a chemical potential for the anyons in combina- 
tion with the second-order term in the first line. There- 



fore, Eq. (50) gets modified as follows: 



AN 



a 



E mj, 

-x,y K ' 



72 t2 £2 
u x' J y u a 



N. 



(55) 



Equation (55) shows that the effect of the off-resonant 
terms becomes negligible if u) a is sufficiently large. How- 
ever, this new contribution to the chemical potential ap- 
pears to lower order in the perturbation expansion than 
the first term (sixth order in the parameters J x ,yi &x,y m_ 
stead of eighth order). Therefore, it is possible to make 
the off-resonant contribution of the same size of the res- 
onant one by slightly relaxing the condition ( 45 1 . Note 
that the second term in Eq. ( 55 ) is always positive and 



can only be canceled if A > 0, which is possible by choos- 
ing suitable frequencies and photon occupations of the 
two modes [see Eq. (48)]. The anyon interaction is there- 



fore repulsive and the final Hamiltonian is similar to the 
case studied in Ref. [THJ 

As a specific example, we plot /i in Fig. [5] as function 
of N, by using in Eq. (55) suitable numerical values of 
the couplings, such that fi is approximately zero. Since 
the dependence on N of Eq. (55) is not exactly linear 
for both terms (due to A and o7 a ), fi cannot be made 
zero identically for arbitrary N. Nevertheless, in a range 
of parameters, it is sufficiently small that the repulsive 
anyon interaction becomes the dominant effect. This we 
show by plotting in Fig. [5] the mean-field gap /i m f, ob- 
tained from Eq. ( 40 ) , which can be much larger than the 



noninteracting gap /!. 
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FIG. 5. Mean- field gap /i m f, given by Eq. (40 1 with A as 



in Eq. (48 1, plotted for different temperatures (solid curves). 
The dashed curve is the bare gap [i of Eq. (55 1. We used 
J z = 1011 K, J X12 = J yi2 = 22 K, (ala x ) = 100, (aja,,) = 0, 
us. 



S x =S y = 0.21 K 



= 3991 K, and ui y = 3980 K. 



The effect of the interactions is very beneficial since it 
allows to obtain /x m f in the mK range, while the origi- 
nal gap without cavities is J x Jy/8J z ~ 30 fiK. However, 
although the interacting regime presents interesting as- 
pects, it seems more useful in practice to take advan- 
tage of the resonant enhancement of the gap discussed in 
Sec. [V] This is realized if each mode couples to all the 
corresponding (x- or y-) links, instead of half of them. 
The final result is the same of Eq. ( 37 1 , including two 
resonant modes 



! 72 

J OL OL 



uj a y 



(aLa a ). 



(56) 



With the same parameters of Fig. [5j Eq. ( 56 ) gives /x ~ 
0.1 K, much larger than the maximum value ^t m f ~ 2 
mK of Fig. [5] As discussed, the effect of the interaction 
becomes larger with system size, but unfortunately the 
perturbative treatment is only strictly applicable within 
a limited range of values of N. Conditions for N are 
discussed in Appendixes |B 2| and[C| and lead to an upper 
bound for N similar to the largest value in Fig. [5j 



VII. CAVITY MODES OUT OF RESONANCE 

In the previous sections, we generally assumed that 
the different cavity modes are resonantly coupled to the 
memory, i.e., uj a ~ 4J Z . We now consider the case of 
largely detuned modes. Since cavity modes with large 
frequency ui a 3> J z have a small effect on the spin model, 
we consider the opposite limit: 



(57) 



and show that many results of the previous sections have 
an analog in this regime. As before, we can still compute 
the SW effective Hamiltonian given by Eqs. (A24) and 



( A26 1 , but it is not possible now to pick a few resonant 
terms at each order of perturbation theory. However, as 
shown in Appendix |DJ the final result can be obtained 
from the calculation without cavity modes, by substitut- 
ing J a ,£c, Ja.y with the corresponding full couplings (J2|: 



Ja 



K„ 



(58) 



Given this simple prescription, the expression Eq. ( 20 ) 



for the Kitaev model with four distinct couplings J Xl2 , 
J yi 2 proves to be useful for several configurations. 

We consider first a single cavity mode coupled homo- 
geneously to the x links as in Eq. ( 35 ) , which amounts 
to set 



■^ci 2 ^ Jx ~t~ &x ~t~ ^ ) ) Jyi,2 ^ 



(59) 



in Eq. (20). As in Sec.|Vj we can neglect the off-diagonal 
terms in the photon Hamiltonian and obtain 



H e s = 



2J 2 



a' a 



16Jf 2^ 

z a 

( 60 ) 

which shows an enhancement of the gap from the cavity 
mode. Even if the gap from Eq. ( 60 ) does not have a 



resonant denominator, it is still possible to increase it by 
driving the cavity to a large photon population. From 
the point of view of the perturbative treatment, the only 
condition to be satisfied in Eq. (60) for the number of 
photons is S x y/ (at a) -C J z - 

We now consider the anyon interactions and choose the 
same coupling scheme of Sec. |VI| with two cavity modes 
a x . y respectively coupled to half of the x- or y links. More 



explicitly, we set in Eq. (|20j) 

J ai -> J a + S a (a a + 



(61) 



where a = x,y. By keeping only terms which are relevant 
for the star and plaquette operators (see Appendix [T5| , 
we obtain 



N 
4X 



Jx^y 



16J| ^ 

z a 
J x Jy $x $\ 
16J3 



/2fl 



a—x,y 



" ( a x al 



16J3 



^ J a S a (a a + oj,) 

+ 4) 



(62) 



which is in complete analogy with Eq. ( |54[ ), except that 
the first ^ a Wa term in the second line (present even 
without cavity modes) is not neglected here, since all the 
energy denominators are non-resonant. In Eq. (62) we 
defined 



N 



(63) 



where the second term is generally a small correction in 
the perturbative regime. As in the previous section, we 
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perform a second SW transformation in the photon field 
and obtain the following effective Hamiltonian by making 
use of the resonant condition cj" ~ u'' 



where 



A = 2 



H 



8J? 



E^ 



uj" 



a',', 



(64) 



(65) 



Similar to Sec. VI we can rewrite Eq. (64 1 in terms of the 



anyon numbers 7i a , in the form of Eq. (39), with A given 
by Eq. (65) above and a chemical potential /i = —AN/2, 



approximately proportional to the system size. We recall 
again that A can be made either positive or negative, 
depending on the sign of {a\.a x ) — (ala y ) and (w" — ui"), 



with A < leading to a stable memory. As in Sec. [Vl 
if we stay a little bit away from the u" ~ u)y resonance, 
the lower-order terms become relevant and we can write 
the following expression for the chemical potential of the 
anyons: 



AN 

'~2 



j2 72 
J x u y 

8J? 



E 



WIN 
16J?w" ' 



(66) 



Secondly, we studied cavity modes coupled homoge- 
neously across the whole memory. We have considered 
both setups with one and two cavity modes, resonantly 
coupled or with very small frequencies. For the case of 
a single cavity we have obtained the dependence of the 
anyon gap from the photon number and frequency of the 
mode, and shown that a large enhancement is possible 
with an excited cavity at resonance. Having two cavity 
modes allows to realize long-range interactions with tun- 
able sign among anyons and a size-dependent chemical 
potential. In this case, the relevant excitation energy is 
a mean-field gap^S which, in some parameter regimes, is 
essentially due to the anyon interactions and is generally 
much larger than the original gap without cavities. Since 
the lifetime of the memory depends exponentially on the 
value of the gap at low temperature, these results might 
lead to a dramatic prolongation of the memory lifetime. 
From a more general perspective, the extension of the 
honeycomb Hamiltonian considered in this work allows 
to have additional control on the properties of the model, 
through the specific features of the coupling scheme and 
the possibility of tuning the parameters characterizing 
the cavity modes. 
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VIII. CONCLUSION 

In this paper we have analyzed a generalization of the 
Kitaev honeycomb model, which includes the influence 
of quantized cavity modes. The coupling is realized by a 
linear modification of the Ising interacti ons w hich should 
be relevant for a variety of systems JM t 2 7 | 30| t38] -yy e h ave 
studied this model in the perturbative regime, when the 
low-energy Hamiltonian is a generalization of the toric 
code and the setup can thus be seen as a quantum mem- 
ory. We have considered different coupling schemes and 
shown that the presence of cavity modes can be useful in 
several ways. 

More specifically, cavity modes allow to read out the 
state of the star and plaquette operators. An efficient 
read-out scheme of the error syndrome is an essential pre- 
requisite to perform any error-correction algorithm I2M5I 
The simplest way to achieve this goal is to couple one cav- 
ity mode to the desired star or plaquette and to detect 
the corresponding frequency shift, similarly to a read- 
out scheme for single su percon ducting qubit in circuit 
quantum electro dynamics J2SI2Z1 A n array of multiple res- 
onators allows to measure all the star or plaquette oper- 
ators simultaneously. 



Appendix A: SW TRANSFORMATION 

For the convenience of the reader, we summa- 
rize here the form alism of the Schrieffer-Wolff (SW) 
transformational 2 !! and show how it can be applied to 
the general class of Hamiltonians of interest in this work. 

1. General formalism 

Consider a Hamiltonian H with a projector Pg on an 
invariant subspace (i.e., a direct sum of eigenspaces) and 
two known projectors P and Q = 1 — P, such that the 
subspace of P has the same dimension of the one of Pjj- 
A SW transformation is defined as a unitary transforma- 
tion U — e s with block off-diagonal S which maps the 
subspace of Ph into the one of P. Therefore, P = e s P H 
and the transformed Hamiltonian H d = UHW is block 
diagonal: 

PSP = QSQ = 0, (Al) 

and 

PH d Q = QH d P = 0. (A2) 
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It is useful to define the superoperator L as 

LA=[S,A], (A3) 

such that, in the superoperator language, the trans- 
formed Hamiltonian H d takes the following compact 
form: 



Ha = e L H. 



(A4) 



Let us now consider an Hamiltonian that can be de- 
composed in a dominant part H and a small perturba- 
tion V = V d + Vodt where the spectrum of H is divided 
in a low-energy space (P projects onto the low-energy 
space) and a high-energy space (Q = 1 — P projects onto 
the high-energy space). By definition, the diagonal per- 
turbation V d and the off-diagonal perturbation V ad sat- 
isfy the following equations 



PV d Q = QV d P = 0, 
PV oA P = QV od Q = 0. 



(A5) 



Since V is assumed to be small, we can expand 5 and L 
in a series, namely 



S =y j S n , and £ = E L n , 



(A6) 



n=l 



where L n A := [S n , A]. If we require that the off-diagonal 
part vanishes at each order, we find the following recur- 
sive relations for the 5„ operators^ (n — 1, 2, 3, ...): 



51 — L Q 1 V od , 

52 = L 1 LiVd) 

1 



L 1 L\H + L L L 2 V d , 



(A7) 
(A8) 

(A9) 



where LqA := [Ho, A]. It is now possible to find general 
expressions for the low-energy effective Hamiltonian at 
each order, namely 



H% = PV d P, 



H is = ^PL 3 V od P - -PLtV od P, (A13) 



H<g = ~PLaV, ld P. 



(A10) 
(All) 

(A12) 



2. Application to a general class of Hamiltonians 

In this section, we want to apply the formalism pre- 
sented above to a general class of Hamiltonians relevant 
to this work: 



H = H a + J2TZ, 



(A14) 



defined on a product Hilbert space %i ® %2- From 
Eq. (17), it is clear that H of Eq. ([I]) has the above 



form. The unperturbed spectrum (given by H ) can be 
described in terms of a single type of excitations in Hi 
(in our case, hard-core bosons) and a set of excitations in 
H2, labeled by an index a (in our case, photons belong- 
ing to different cavity modes). The operators Tjj acts on 
Hi by changing the number of hard-core bosons by n (in 
our case, n = 0, ±2) and on Hi according to the index j. 
This can assume the values j = ma, indicating that the 
number of photons of type a changes by m (in our case, 
m = 0, ±1). As discussed in the main text, we also use 
the convention 0a = 0, such that there is a unique set of 
T° operators (which do not change the photon number 
of any mode). 

At this point, we define the low-energy space of this 
system as the lowest energy eigenspace in Hi- Therefore, 
P is the projector onto the subspace with no hard-core 
bosons, and Q = 1 — P onto the subspace where some 
hard-core boson is present. The diagonal and off-diagonal 
perturbations are given by 



(A15) 



= :B 



V od = PBQ + QBP (A16) 
We assume in the following, as appropriate for our case, 



PTJ = %P = 0, 



(A17) 



and that the energy change due to the T£ operators is 
given by 5e J n . Therefore, the following commutation re- 
lations hold: 



[H ,Tl]=6elTl 



(A18) 



In the simplest case, the excitations also correspond to 
energy quanta, as in Eq. (18). 



It is now possible to calculate 51,2,3 according to 
Eqs. ( |A7]|A9l ), which give 



-H.c. 



5 2 



5, 



EE 

i,j n,m 

EE 



PT l n Q 
PTtQTjQ 



H.c. 



(A19) 
(A20) 



PTiQTlQTfQ 



1 PTlQTfPTiQ + PT^QTiPTfQ 

2 PT^QT^PTfQ 



Se?) 



H.c. 



(A21) 



j.ri 



where every term with n — vanishes because of 
Eq. (A17) and H.c. denotes Hermitian conjugates (note 
that all the S p are anti-Hermitian) . These Hermitian 
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conjugate terms can be written in a form similar to the 
ones explicitly appearing in the square parenthesis, by 
making use of 



(Tn) 



and 5e_ 



-del 



(A22) 



With the help of Eqs. ( |A10|A13[ ), we finally obtain the 
low-energy effective Hamiltonian up to fourth order: 



H 



H. 



(i) 

eS 

(2) 
off 



JEE 



i,j n,m 



H.c. 



(A23) 
(A24) 



H. 



(3) 

eg 



PT*QTiQT?P 



rr(4) _ 
^eff - 



i,j,k,r n,m,l : q 



H.c. 



SeiiSei + Sein) 

2 -PTlQTlPTfQTZP 



(A25) 



1 PTiQT l k PTlQT^P + PTj n QT^ PJf QT r q P 
1 PT^QT^PTfQTqP + ZPT l n QT r q PT^QT^P 
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Se^SeLSef 



+- 



PTn QT] n QT t k QT q P 



H.c. 



(A26) 



where again H.c. are Hermitian conjugate terms. 



Appendix B: SW TRANSFORMATION WITH A 
SINGLE RESONANT CAVITY MODE 



We consider in this appendix the application of the 
SW formalism to the case when a single cavity mode is 
present, coupled only to x links. This is useful for both 
Sees. IV and |Vj where a read-out scheme and the effect 
of the cavity mode on the gap were studied. 



1. Derivation of the effective Hamiltonian 

In this case, the perturbation is given by T° and 
operators, where we drop the index a since only one mode 



is present. The treatment of the previous Appendix A 2 
can be easily applied to this case, where 5e J n of Eq. ( A18 ) 



is as in Eq. ( |18[ ), with uj a = uj. Since we are interested in 
the resonant case (4J Z — uj) <C 4 J z , uj, the relevant energy 
differences are 



5e 



+ 2 



-8J 



4J Z 



6e° = 0. 



(Bl) 
(B2) 



are obtained from Se +2 , <5el 2 , an d <5cq, which leads to 
1 fl 



H c ff = uja'a ■ 



+ 



1 



{AJ z -ujf V« 



4 J z - uj V 2 
•5 



T+ 2 T +2 + H.c. 



1, 



r p-r r p— r Y ] ~r T 1 n^' T 7 

- i -2 J +2 i -2 i +2 - 7 J -2 J -2 J +2 J +2 



■V 2 T +2 T± 2 T+ 2 -1t+T T °T +2 + H.c. ) . (I!:!) 



where T +2 



77 2 + T"^ as defined in Eq. \17\ 



The effective Hamiltonian (B3) is valid in the subspace 
with zero hard-core bosons, but is not diagonal in the 
number of photons. Therefore, it represents an effective 
Hamiltonian for the cavity and it requires to be further 
diagonalized to obtain the energy eigenstates. The ex- 
plicit form of the photon Hamiltonian depends on the 
specific coupling scheme of the cavity to the links of the 
model, which is reflected in the operators. As an 
example, we assume the homogeneous coupling scheme 
(35), to all the x links of the model, and consider first 



the terms up to second order. Evaluating the first line of 
Eq. (B3| in terms of effective spins and photon operators 



gives 



H, 



eff 



SIN 



-a' a - 



(a + at) ^ 
4 J z - uj v ' 2 



,t2* 



N 



4J 2 



(B4) 



The first two terms appear in Eq. (36) while the off- 



diagonal terms can be eliminated with a second SW 
transformation, in the photon operators only. Such off- 
diagonal terms involve excitations with energy uj or 2uj 
and therefore give small corrections in the resonant limit 
we are interested. Therefore, we have neglected them in 
Eq. (l36|). 



Evaluating all the fourth-order terms appearing in 
Eq. (B3 1 is cumbersome, but this is not necessary for our 
purposes. As explained above, we keep only the terms 
already diagonal in the photon number, while the other 
ones give only small nonresonant corrections. Therefore, 
we can simply substitute all the T n operators in Eq. ( B3 ) 
by T~ . Furthermore, we are interested to the coupling 
of the photon mode to the W a operators, which are the 
product of spin operators involving four distinct links of 
the model (of which two are x links and two y links). To 
understand the relevant terms at fourth order we note 



that, for the specific coupling scheme (35), only the x 



links are coupled to the cavity mode. Therefore: (i) the 
T° operators involve a sum of terms with all the links of 
the model, while (ii) the operators are sums of terms 
with only x links. Since only the Tq operators contain 
terms relative to the y links, it becomes clear that only 
the T_ 2 T§TqT+ 2 com bination can couple to the W a op- 
erators. Finally, (B3) can be simplified to the following 
form 



It is then possible to directly apply Eqs. ( A24) and ( A26 1, 
keeping only the terms where the energy denominators 
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annihilates the two hard-core bosons. Therefore the links 
involved form a single closed loop, of the form abed (if 
all four links are different, which contributes with a term 
proportional to W a ) or abba (if the links are repeated, 
which contributes with a term without W a ), see Fig. [6] 

On the other hand, the remaining three fourth-order 
terms appearing in Eq. (B3), with their Hermitian conju- 



FIG. 6. Pictorial representation of the connected and discon- 
nected terms entering the perturbation expansion ( B3 1 . Both 



loops abed and abba represent the different contributions aris- 
ing from the connected term T^ 2 TqTqT^ 2 . We see that only 
abed loops lead to terms proportional to since all four 
links are different. This is not the case for the abba loop 
and the corresponding terms are constants. The remaining 
fourth-order terms entering Eq. ( B3 1 can be represented by 
two disconnected loops aa and bb. Since the number of dif- 
ferent positions of each aa and bb loop is proportional to TV, 
the contribution of the disconnected terms is ~ TV 2 . 



gates, are disconnected. To see this, we note that at most 
one of the four operators, the T± 2 , can act on a y-link, 
while the other three, of the form , only contain x-link 
terms. Therefore, the only possibility is that two pairs 
of hard-core bosons are created on distinct x links and 
annihilated on the same two links, a processes that can 
be represented with two small disconnected loops aa and 
bb, see Fig. [6] These processes can contribute with terms 
which are of order TV 2 instead of TV. Clearly, this does 
not occur for the read-out schemes discussed in Sec. [TV 



where each cavity mode is coupled locally to a single 
plaquette or star, but is relevant for Sec. |Vj By assuming 
the homogeneous coupling scheme (35 ), the disconnected 
terms of Eq. (B3) evaluate to 



which evaluates to Eq. (36), by expressing it in terms 
of photon and effective spin operators. As it turns out, 
the same reasonings apply to the read-out methods we 
discussed in Sec. |IV| where as well only couplings to the 
x links appear. Therefore, Eq. (B5) is also the relevant 



\ S (ij - Jxa ^ a + ^ at2fl2 + + H.c.) 

5 5lN 2 



(J x a'aa' + 5 x a) aa) a + 5 x a)aa' 2 + H.c.) 

( J x a' 2 a + 5 x a) aa) a + S x a' 3 a + H.c). 

(B6) 

@, which evaluate to Eqs. ([27} and @ respectively. It is immediate to check that, if the photon operators are 



8 (4J Z - wi) 3 
1 S X N 2 
8 (4Jz - Wl ) 3 



effective Hamiltonian for the coupling schemes (26) and 



Note that we have always neglected in the final results 
a subleading constant correction of the cavity frequency, 
i.e., the a^a term appearing in fourth order and indepen- 
dent of the Wa operators. 



2. Disconnected terms 

We would like to comment here on a special aspect 
which arises in the presence of cavity modes, i.e., the ap- 
pearance in the perturbative expansion of non-extensive 
terms in the size of the system TV. These originate 
from disconnected terms in the perturbation expansion 
which, if the Hamiltonian is local, give a vanishing 
contribution.^ This result, known as linked-cluster the- 
orem, is not valid here because of the long-range nature 
of the coupling of the photon modes, which are delocal- 
ized and can interact simultaneously with links at any 
distance across the lattice. 

For the case of interest here, the only connected fourth- 



(t) 



(*) 



Eq. (B6) 



replaced by complex numbers, i.e., a 
gives an extensive contribution since all the N 2 terms 
cancel, in agreement with the linked-cluster theorem. If 
we use the bosonic commutation relation [a, a r ] = 1, 



then we can simplify Eq. (B6| to the following form 



3TV 
"~2~ 
N 2 



^[5 x (a t2 a 2 + a t2 aa f + H.c.) + J x {a) 2 a + a f a 2 )] 



order term appearing in Eq. ( B3 1 is the last one, with 
its Hermitian conjugate. Consider in fact the action of 
T+ 2 T§TgT +2 : the first operator from right to left, T +2 , 
creates two hard-core bosons on a certain link. The sec- 
ond and third operators make such hard-core bosons hop, 
and therefore they have to act on links where only one 
hard-core boson is present. Finally, the last operator T_t 2 



+ — 5 x [86 x a^a + 5J x (a + a 1 ) + 6S x {a 2 + a f2 )], (B7) 
8 

where the TV 2 contribution (second line) does not vanish. 

This result does not represent a fundamental problem 
since, for a finite system, terms proportional to TV 2 or 
higher powers of TV can certainly exist. However, the 
final perturbative expressions are not extensive quantities 
in the thermodynamic limit TV — > oo. This behavior can 
be attributed to the failure of the perturbative treatment 
since, besides the conditions J x , y , 8 x yJ (at a) -C (4J Z — cj), 
the requirement to have a power series in terms of small 
expansion parameters also leads to an upper bound on 
the size of TV, namely 



(B8) 



5 X VN«. {U z -u) x ). 



This condition is not very restrictive, since the cou- 
pling S x is expected to be small. Furthermore, an in- 
teresting limit is obtained for a vanishing coupling S x 
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and a large number of photons (a^d). In this case, for 
system size N <C (a^a), the non-extensive terms can 
be generally neglected in comparison to the extensive 
ones. For example, in Eq. (B7| this condition would give 
S^N(a 1 'a} 2 S> SlN 2 (a^a), by comparing the first and sec- 
ond line. 



Appendix C: SW TRANSFORMATION IN THE 
PRESENCE OF FOUR DIFFERENT LINKS AND 
TWO RESONANT CAVITY MODES 

In this section we present how the SW transformation 
is applied to the scheme of Sec.|Vl]for which the four links 
belonging to each W a operator have all distinct couplings, 
see Eqs. @, {42]), and Fig. g By using Eqs. ( |A24| and 
(A26) one can obtain the zero hard-core boson effective 
Hamiltonian up to fourth order. Since we are interested 
in a regime where (4J Z — u> x ,y) ^ ^x,yi we keep the res- 



onant terms only. Furthermore, it is sufficient to keep at 
fourth order only the terms that couple to the star and 
plaquette operators. This leads to: 



H, 



eff 



E 



T +a T T 0( T 



T« 2 (T° 2 +T« 2 +T- 2 a )+H.c. 



+2 



• T 



+ 2 



2(4J 2 
f H.c. 



U) a ) 



2(4J Z 



where a = y if a = x and vice versa. In the above ex- 
pression, terms like T^ 2 T^ 2 Tl 2 Tl 2 or T^T^T^T^ do 
not appear, because they do not contribute to the reso- 
nant coupling to the W a operators. It is worth mention- 
ing that the conditions in order to treat the Hamiltonian 
within perturbation theory are the same as in Appendix 
|B} namely, 



S a y (ata a ), S a VN -C min(4J 



,)• (C2) 



By writing the Hamiltonian (CI I explicitly in terms of 



effective spin operators and photon operators, we obtain 
the following zero hard-core boson effective Hamiltonian 



H, 



eff= 



E 

a—x,y 

si 



ui' a a\a a 



4J Z - cu a 
N 



4) 



N 



4J Z — uj a 
J 2 J S 



a) 2 )- 



(4J Z 



4)E^ 



(C3) 



In the strictly resonant case [see Eq. (45 )] we can drop all 
the nondiagonal terms except the last line, which leads 
to Eq. (44). If we want to calculate the leading nonres- 



terms in the first and third line, 



keep also the [a a + a) a 

which leads to Eq. |54| ). This is because, as explained 
in Sec. |VI C\ these two nondiagonal terms combine to a 
higher-order ^ a W a contribution, diagonal in the pho- 
ton modes, after a second SW transformation. On the 
other hand, (a 2 a 4- a* 2 ) does not appear together with 
S a anc ^ therefore can be dropped, for the purpose of 
obtaining corrections to the star and plaquette couplings. 
Finally we note that the third line, after the second SW 
transformation in the photon operators, also leads to a 
interaction term (J^ a W a ) 2 , but this is generally negligi- 
ble compared to the one from the last line of Eq. ( C3 ) 
and is not considered in Sec. IVI CI 

As a last note, a conservative requirement to perform 
the SW on Eq. ( C3 ) is that that the off-diagonal opera- 



tors are much smaller than the corresponding gaps. From 
the first line of Eq. ( C3 ) we obtain 



J n S, 



4 J z - u a 
and from the fourth line 

JxJ y SxSy\J (a x ai)(a y ay) 



2(4J Z 



L0 o 



N^\u' x -uj'\. 



(C4) 



(C5) 



(CI) Finally, the condition from the third line of Eq. (C3| is 



less restrictive than Eq. (C4), because of an additional 



small factor J|/(4J Z — uj a )' z . These conditions, together 



with Eq. \C2\ , pose a restriction onto the size N at which 
the perturbative treatment is justified. 



Appendix D: SW TRANSFORMATION IN THE 
PRESENCE OF SMALL FREQUENCY MODES 

The main purpose of this section is to prove the pre- 
scription (58) of Sec. VII to obtain the effective Hamil- 



tonian in the presence of off-resonant cavity modes. 

First we note that, if no cavity is present, only the T° 
operators are nonvanishing and T n = T®. Therefore, we 
can drop in Eqs. ( A23||A26 ) the summations in the pho- 
ton indexes i, j, k, r and set all the operators equal to 
T n . On the other hand, by making use of Eq. (57) for 



the case with small frequency modes, we can approximate 
5e l n ~ 2nJ z , i.e., neglect the photon frequency shifts in 
the denominators of Eqs. ( A23f A26 ). This allows to per- 
form the summations in the photon indexes and express 
Eqs. ( A23||A26 1 in terms of the T n operators, as for the 
case without photon modes discussed above. The dif- 
ference is that the T n are now photon operators, simply 
obtained from the ones without cavities by substituting 
Ja.,k ^ > -Ka.fc (where k — x,y). Furthermore, it is also 
clear from Eq. S that the K a k are all commuting oper- 
ators, and therefore they can be treated in the same way 
of the J a ,fe in the derivation of the effective Hamiltonian. 
This shows that it is sufficient to apply Eq. (58), if the 



onant correction to the chemical potential, we have to effective Hamiltonian without cavity modes is known. 
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By applying this prescription to the second coupling 
see Eq. (61 )], we obtain the following 



scheme of Sec. 
effective Hami 



VII 



toman: 



4./ 



'q'- 
16J? 



(a a +al)N 

a—x,y 

(a a +< )— - 3 (aa + aL)2^ 



2 t2 



where ui'^ are defined in Eq. (63). In Eq. (Dl|, we only 



kept terms from the first and last lines of Eq. ( 20 ) while 



dropping all the fourth order terms which do not couple 
to the Wa operators and some inessential constants. 



We note now that the {a 2 a 



rt2 



) term, in the second 
line of Eq. (Dl|, does not have a corresponding £^ a W a 
term, and that the off-resonant coupling to the W a op- 
erators in the last line does not have a corresponding 

{a x a y + a\.a!y) term at lower ordei\ Therefore, after a 

second SW transformation of Eq. (Dl ), these two terms 



give negligible corrections to the effective Hamiltonian 
for the star and plaquette operators and we have dropped 



(Dl) them in Eq. (62) 
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